Binary Pulsar Constraints on the Parameterized post-Einsteinian Framework 
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We constrain the parameterized post-Einsteinian framework with binary pulsar observations of 
orbital period decay due to gravitational wave emission. This framework proposes to enhance the 
amplitude and phase of gravitational waveform templates through post-Einsteinian parameters to 
search for generic deviations from General Relativity in gravitational wave data. Such enhancements 
interpolate between General Relativity and alternative theory predictions, but their magnitude must 
be such as to satisfy all current experiments and observations. The data that currently constrains 
the parameterized post-Einsteinian framework the most is the orbital period decay of binary pulsars. 
We use such observations to place upper limits on the magnitude of post-Einsteinian parameters, 
which will be critical when gravitational waves are detected and this framework is implemented. 



Introduction. Gravitational waves (GWs) will allow us 
to learn about the gravitational interaction in regimes 
that are currently inaccessible by more conventional, 
electromagnetic means. Binary black hole and neutron 
star mergers, for example, lead to gravitational fields 
that are intensely strong and highly dynamical, a regime 
where General Relativity (GR) has not yet been tested. 
GW theorists and data analysts will need to be able to 
make quantitative statements about the confidence that 
a certain event is not just a GW detection but one con- 
sistent with GR. 

The parameterized post-Einsteinian (ppE) frame- 
work [l| was devised precisely for this purpose: to search 
for statistically significant GR deviations or anomalies in 
GW data and, in their absence, to quantify the degree 
of belief that a GW event is purely described by GR. 
This framework enhances the waveform templates used 
in matched-filtering through parameters that character- 
ize GR deformations. In practice, this is achieved by 
adding to GR templates amplitude and phase correction, 
with magnitudes depending on certain ppE parameters. 

Any framework that modifies GR must comply with 
Solar System and binary pulsar observations. These mea- 
surements already strongly constrain GR deviations in 
weak and moderately strong fields. The ppE framework 
was constructed on a maxim of compliance with current 
observations, which can be enforced by requiring that 
the magnitude of the ppE correction be such as to sat- 
isfy current constraints. Until now, this maxim had not 
been quantitatively enforced because it was thought that 
it would be difficult to relate the ppE deformations to 
Solar System or binary pulsar observations. 

We have here found a relatively simple way to relate 
the ppE framework to current experiments. As shown 
in [1|, modifications to the dissipative and conservative 
sectors of the theory lead structurally to similar ppE cor- 
rections to the waveform. We find that to constrain the 
ppE framework with current experiments at least ini- 
tially, it suffices to consider dissipative corrections only, 
while keeping the conservative sector unmodified. Such 
dissipative corrections modify the amount of orbital bind- 
ing energy carried away by GWs, which affects directly 



the orbital period decay in binary pulsars. 

The relatively recent discovery of the binary pulsar, 
PSR J0737-3039 0, has provided particularly powerful 
GR tests [3]- This pulsar is highly relativistic, with an 
orbital period of about 2 hours, and has an orbital geome- 
try favorable for measuring quantities such as the Shapiro 
delay with sub-percent precision. Such data has been re- 
cently used to constrain alternative theories of gravity to 
new levels [J|. 

In this paper, we relate such sub-percent accurate mea- 
surements of the orbital decay of PSR J0737-3039 to 
constrain the ppE framework and its templates. Be- 
cause of the structure of the ppE correction to GWs, 
these constraints are relational, i.e., they are of the form 
ItI/"^ < F{5, A), where (7, c) are ppE parameters, / is the 
GW frequency and F{5, A) is some function of the accu- 
racy 5 to which the orbital decay has been measured and 
system parameters A, such as the mass ratio and total 
mass of the binary. Thus, given a value for c, the mag- 
nitude of 7 is constrained by binary pulsar observations 
to be less than some number related to 5, f and A. The 
relational constraint found in this paper will be crucial 
in the implementation of the ppE framework in a realis- 
tic data analysis pipeline once GWs are detected. In the 
rest of this paper, we follow mostly the conventions of Q 
with geometric units G = c = 1. 

Basics of the ppE Framework. The main GW observ- 
able is the so-called response function, which describes 
how an interferometer reacts to an impinging GW. In 
GR, this function is given by 



h^^{t)=F+h%^{t) + F^h^,^{t). 



(1) 



where i^+.x are beam-pattern functions and h'^\ arc the 
plus and cross GW polarizations, built from contractions 
of the metric perturbation with certain polarization ten- 
sors 0. For quasi-circular binaries, these polarizations 
can be Fourier transformed in the stationary-phase ap- 



proximation [6|-|8| to yield 

M u2/3 



= (1 + COS^ i) e-'(*OR+V4-2/3) 



where (t, /3) are the inchnation and polarization angles, 
Dl is the luminosity distance from source to observer and 
F is the rate of change of the orbital frequency due to 
GW emission. This frequency is defined as F = (l/27r)\i', 
where \1/ is the orbital phase, and it is also equal to half 
the Fourier or GW frequency /, i.e., F = f/2. The 
quantity ^Pgr is the GR GW phase in the Fourier domain, 
which can be computed via 



^oM) = 2nf'C(2-l,^dF'. 



(3) 



The quantity u = nAif is a dimensionless frequency pa- 
rameter, where A4 — rf'^m is the chirp mass, with 77 — 
17111712/171^ the symmetric mass ratio and m = mi + TO2 
the total mass. From Eq. ([T]), it follows that the Fourier 
transform of the response function in the stationary 
phase approximation is simply /iqr = F^h^ + Fxh°^. 

The ppE framework proposes that one enhances the 
GR response function via an amplitude and a phase cor- 
rection. In the Fourier domain and in the stationary 
phase approximation, one can parameterize the response 
function for a GW from an unequal-mass, binary, quasi- 
circular inspiral as [l|, |9| 



h^^{l + a7^-u-)e^^^"-\ 



(4) 



where (a, c, a) are ppE amplitude parameters and (/3, d, b) 
are ppE phase parameter. Such a correction arises gener- 
ically if one modifies F — E {dE\^/dF)~^, which in turn 
can arise either due to a modification to the GW lumi- 
nosity E (the dissipative sector) or to the orbital binding 
energy £'b (the conservative sector). As explained in [l|, 
this degeneracy breaks the one-to-one mapping from a 
ppE waveform modification to a specific alternative the- 
ory, as one cannot tell whether the change arose in the 
dissipative or conservative sector. 

Gravitational Wave Luminosity . We now compute the 
energy carried by ppE GWs. As is clear from Eq. ([2]), 
the GW amplitude depends on F, which by the chain 
rule can be related to E as explained below Eq. ^. We 
can construct E directly from /i+ or /ix via 



E - -f -Dl/gr 



dn 



h. 



(5) 



where /gr is the rate of change of the GW frequency, and 
dQ — sin t di dp integrates over the (i, /3) dependence of 
the waveform. Notice that Eq. ^ agrees with Eq. (2.38) 
in @. Substituting for h using Eq. ([5]), we find 



E = Ea^\l + aifu'' 



(6) 



where -Egr is the GR expectation for the GW luminosity: 



^GR = -^fonfDl I dn 



hi"" 



K 



(7) 



One can also obtain an expression for the GW luminos- 
ity in terms of the GW phase only, as this also depends 
on F as shown in Eq. ([3]). Noting that d'^'^/dp = 7ri^"\ 
we can write the GW luminosity as 



d^^i> 



Since * = vj/^^, + /Jry'^u'', we find that 



(8) 



E = Ea 



l + Ti^M'^l3r]'^h{h-l)u 



b-2 



dp 



(9) 

where Sgr can be written in terms of the GW phase as 

(10) 



^GR--g/GR-M U ~Jf2- 



This shows that measurement of E from a circular bi- 
nary pulsar would allow us to constrain both the ampli- 
tude and the phase ppE parameters. All known binary 
pulsars, however, are in eccentric orbits. The /qr used in 
Eqs. ([5]) and (|8]), or equivalently the i^cR used in Eqs. (E]) 
and ^, must be that of an eccentric orbit, namely [lOi] 



Ea 
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24 



96 



(11) 

The ppE corrections to E computed in Eqs. ^ and (|9]) 
are built from circular, non-spinning ppE templates, the 
only ones currently available. The ppE corrections to the 
GW luminosity computed here are really only the lead- 
ing order terms in a post-circular expansion [8|, i.e., an 
expansion in small eccentricity e. If the eccentricity is 
smaU, such as with PSR J0737-3039 (e = 0.088), such 
an approximation is well-justified and will not strongly 
affect the constraints we place on the ppE framework. 

Orbital Period Decay. The GW luminosity enters 
into binary pulsar observables through the orbital decay: 
P/P = (3/2) Eb/Eb = -(3/2) E/Eb, where in the second 
equality we used energy balance: the amount of binding 
energy lost by the system is equal to minus the amount 
of energy carried away by GWs, Eb = —E. Using Eq. © 
and dl]), we then find that the P corrected by amplitude 
ppE parameters is 



\ / GR 

while that corrected by phase ppE parameters is 



(12) 




A + 15 /3 ^rf 5(6 _ 1) u'^+'^/A . (13) 



The quantity (P/P) gr stands for the orbital decay in GR 
for an eccentric inspiral, namely [10] 



96 77771 



a4 ^ ^ V 24 96 



(14) 

Recall again that the ppE corrections [the second terms 
inside the parenthesis of Eqs. (IT^ and P^ ]. are only 
valid to leading order in the post-circular approxima- 
tion. In deriving these expressions, we have used the fact 
that the observed P/P is very close to the GR value: 
(P/P)obs = (P/P)gr(1 + 5)- The observational error 
5 = {5P)/P <^ 1, meaning that the error on P domi- 
nates over the error on P. 

Since binary pulsar observations have confirmed GR 
up to observational error, we can now place relation con- 
straints on the ppE framework. Focusing first on the 
amplitude ppE parameters, we find that 



l«l< 



1 5 



For the phase ppE parameter, 

1/31 < ^ ' 



48|6| \b - 1\ rj^u^+^Z^ 



(15) 



(16) 



A binary pulsar measurement of P to an accuracy 6 
allows us to constrain a and /3, given some value for 
(a,b,c,d), the symmetric mass ratio and the GW fre- 
quency, or equivalently, the orbital period. 

Before proceeding, let us first discuss the apparent de- 
generacy between the amplitude and the phase correc- 
tion. Comparing Eqs. (fT2|) and ([T3l) . one realizes that 
if changes to the GW amplitude and phase are due to 
the same mechanism (for example, if only E is mod- 
ified), then we must have a = b + 5/3, c = d, and 
/3 — 5a/[48b{b — 1)]. The ppE scheme, however, al- 
lows for modifications to both the dissipative (E) sec- 
tor and the conservative (Ei,) sector. If both sectors are 
modified, there will be two sets of independent modifi- 
cations, one to the phase and one to the amplitude. If 
a ppE correction is introduced to the GW amplitude, 
then it is constrained by Eq. ((T^); if a ppE correction 
is introduced to the GW phase, then it is constrained 
by Eq. (fT5)) . These constraints on the amplitude and 
phase ppE parameters are thus independent from each 
other, even though a constraint or measurement of them 
would not allow a one-to-one mapping to a conservative 
or dissipative modification. Thus, conservative and dissi- 
pative modifications are in fact degenerate, even though 
the phase and amplitude measurements are not. 

Binary Pulsar Constraint. Let us now employ the re- 
cent measurements of 2, 3i] on PSR J0737-3039 to con- 
strain (a,P). This binary consists of two neutron stars 
with component masses mi — 1.3381(7) M© and 7712 = 
1.2489(7) M0 in an almost circular orbit with eccentric- 
ity e = 0.0877775(9) and period P = 8834.535000(4) s. 



The symmetric mass ratio is 7/ ~ 0.24970, the chirp mass 
isM~ 5.5399 X 10"*^ s, the GW frequency is / = 2/P ~ 
2.263842976 x 10^"* Hz, and the reduced frequency is u ~ 
3.940046595 x 10~^. The time derivative of the period is 
measured to be P = —1.252(17) x 10^^^, which implies an 
uncertainty of 5 = 0.017 x 10"iV(l-252 x lO^^^) ^ ^q^^ 
This uncertainty is comparable to the systematic error 
in the ppE parameters due to the neglect of eccentric- 
ity effects in PSR J0737-3039, which roughly scale as 
e^ ~ 0.0077. An increase in the accuracy of the P 
measurement, reducing (5, would not allow us to place 
stronger constraints until the ppE templates are extended 
to include eccentricity. 

Figure [T] plots the double binary pulsar constraints on 
(|q;|, |/3|) as a function of the exponent ppE parameters 
(a, 6) for fixed (c, d). The area above the curves is ex- 
cluded by binary pulsar observations, forcing (a, /3) to be 
smaller than a value which depends on (a, b, c, d). Gener- 
ally, if a < —0.4 then \a\ < 10^^ for all plotted values of c, 
while if 6 < -1.9 then |/?| < 10"^ for all plotted values of 
d. For a > 0.2 and b > —4/3, a and /? can be greater than 
unity for all plotted values of (c, d). This makes sense: 
as (a, b) become large and positive, the ppE correction 
becomes smaller for low reduced frequency sources. 

These constraints are consistent with other constraints 
on GR deviations from binary pulsars. For example, 
one can place a generic constraint on the time-variation 
of Newton's constant G with a binary pulsar observa- 
tion [ll|,|I2: G/G < ((5P)/(2P), where 5P is whatever 
part of P that is otherwise unexplained. Using PSR 
J0737-3039 [i,i] one infers that G/G < 3 x 10"" yr-^. 
Allowing for Newton's constant to be a linear function 
of time leads to a modification that can be mapped 
to Eq. gj with |a| = (5/512)(G/G)M, c = 3/5 and 
a — —8/3 for the amplitude parameters, and |/?| = 
(25/65536) (G/G)M,^d = 3/5 and b = -13/3 for the 
phase parameters [1J| . From the binary pulsar constraint 
on G/G, we then infer that |a| < lO'^^ and \/3\ < lO'^^, 
which is consistent with Eqs. P^ and ^TE\i and Fig.[TJ 

Our constraints on a look extremely strong (e.g., for 
a < —2, then \a\ < 10~^°). However, this does not im- 
ply that the unconstrained region (below the curves in 
Fig. [1]) is uninteresting. For example, constraining G/G 
below 10~^^ yr~^ or 10~^^ yr""'^ implies constraining |a| 
below 10~^^. This is interesting as there are GR modi- 
fications that suggest G/G deviations of this order may 
be present J^]. On the other hand, the smallness of the 
y-axis of Fig. [1] does suggest that a and /3 are perhaps 
not the best "coordinates" with which to measure GR 
deviations when a and b are sufficiently negative. 

We conclude this discussion with some caveats on the 
constraints we find. First, we have here neglected the 
effect of eccentricity in the ppE correction, although this 
is accounted for in the GR part of P/P. We have not 
studied eccentric ppE templates because these do not yet 
exist. This is due to the difficulty in constructing analyti- 
cally simple Fourier transforms of eccentric inspiral wave- 
forms in the stationary phase approximation [S,]- Second, 
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FIG. 1: Left: Constraint on |a| as a function of a for fixed c. Right: 
below the curves is allowed, while the area above is ruled out. 



Constraint on |/3| as a function of h for fixed d. The area 



because of the ppE degeneracy between conservative and 
dissipative corrections, we have only examined dissipa- 
tive ones here; recall, however, that amplitude and phase 
measurements are truly independent. Since the conser- 
vative sector can be thought of as unmodified, this allows 
us to use the GR measured values for the components' 
masses, as the Shapiro time delay and periapsis preces- 
sion are the same as in GR (leading to identical results for 
7711 and 7712 in GR and in the ppE extended theory). We 
could have instead allowed for both conservative and dis- 
sipative corrections, and then analyzed how these affect 
all binary pulsar observables. A combined analysis of all 
these effects will presumably lead to a stronger bound on 
the ppE parameters; we leave this to future work. We 
note that we could have studied constraints on ppE from 
Solar System observations. However, the exquisite accu- 
racy of the double binary pulsar measurements, and the 
fact that this is a much stronger-field source than any 
Solar System one, means that Solar System constraints 
will not be as stringent as the ones discussed here. 



Implications for GW Data Analysis. Once GWs are 
detected, one would like to implement the ppE frame- 
work in a realistic data analysis pipeline. Such a pipeline 
will likely employ techniques from Bayesian analysis |15{ , 
which relies heavily on the priors chosen for the param- 
eters searched over. The prior tells us whether certain 
regions of parameter space are allowed or likely to oc- 
cur in Nature. The priors for the ppE parameters should 
be constructed following current Solar System and binary 
pulsar constraints. Equations (|T5]) and (TlGl) represent the 
most stringent prior found to date for these parameters 
using binary pulsar observations. 
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